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ABSTRACT 

Dust is an important tracer of chemical enrichment in primeval galaxies and it has also important implications for their evolution. 
So far, at z >6, close to the re-ionization epoch, the presence of dust has been only firmly established in quasar host galaxies, which 
are rare objects associated with enormous star formation rates. The only non-quasar host galaxy, with modest star formation rate, 
for which dust extinction has been tentatively detected at these early cosmic epochs, is the host of gamma ray burst GRB050904 at 
z = 6.3. However, the claim of dust extinction for this GRB has been debated in the past. We suggest that the discrepant results occur 
primarily because most of previous studies have not simultaneously investigated the X-ray to near-IR spectral energy distribution of 
this GRB. The difficulty with this burst is that the X-ray afterglow is dominated by strong flares at early times and is poorly monitored 
at late times. In addition, the Z band photometry, which is the most sensitive to dust extinction, has been found to be affected by strong 
systematics. In this paper we carefully re-analyze the Swift/XRT afterglow observations of this GRB, using extensive past studies of 
X-ray flare properties when computing the X-ray afterglow flux level and exploiting the recent reanalysis of the optical (UV rest 
frame) data of the same GRB. We extract the X-ray to optical/near-IR afterglow SED for the three epochs where the best spectral 
coverage is available: 0.47, 1.25, and 3.4 days after the trigger. A spectral power-law model has been fitted to the extracted SEDs. 
We discuss that no spectral breaks or chromatic temporal breaks are expected in the epochs of interest. To fit any UV rest-frame dust 
absorption, we tested the Small Magellanic Cloud (SMC) extinction curve, the mean extinction curve (MFC) found for a sample of 
QSO at z > 4 and its corresponding attenuation curve, as well as a starburst attenuation curve, and the extinction curve consistent with 
a supernova dust origin (SN-type). The SMC extinction curve and the SN-type one provide good fit to the data at all epochs, with an 
average amount of dust absorption at A„si = 3000 A of A3000 = 0.25 + 0.07 mag. These results indicate that the primeval galaxy at 
z = 6.3 hosting this GRB has already enriched its ISM with dust. 

Key words, dust, extinction - y-ray burst:individual: GRB 050904 



1. Introduction 

. The presence of dust at high redshifts (z>6) is fundamental both 
' for the formation and evolution of the stellar populations in early 
" galaxies, as well as for their observability. High dust masses 
at z>6 have been detected in the host galaxy of some power- 
ful quasars, through de tection of intense mm/submm (far-IR rest 
. frame) d ust emission (iBertoldi et al.1 120031: iBeelen et ani2006t 
■ I Wang et al. 2008). However, these powerful quasars are very 
] rare objects, probably lying in extreme over-density peaks in 
. the early universe, and hosted in extremely powerful startbursts, 
then forming stars at a rate higher than 1000 M^/yr. It is impor- 
tant to understand whether dust is also present in less extreme 
galaxies at z>6, which are more representative of the primeval 
' galaxy population. 

The main issue in the formation of dust in the early uni- 
verse is the time scale for dust production mechanisms. In the 
local universe dust is mostly produced in the atmosphere of AGB 
stars, which require timescales ranging from a few 100 Myr to 
Gyrs to enrich the ISM with significant amounts of dust. These 
timescales may prevent young primeval galaxies in the early uni- 
verse to have significant amounts of dust (but see also Valiante et 
al. 2009). A lternatively, dust prod uction can occur in the ejecta 
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of SNe (e.g.' Todini & Ferrarall200 1). which can enrich the inter- 
stellar matter (ISM) on a very sho rt timescale, but which can 
also destroy dust through shocks (iBianchi & Schneider! 120071: 
iNozawa et aDl2007l) . 



If the dominant dust production mechanism changes over 
the cosmic epochs, this may result in a change in the dust 
properties (e.g. traced by the extinction curve) as a function 
of redshift. In the local universe and at low redshift, differ- 
ent systems are characterized by different extinction and at- 
tenuation curves. The SMC extinction curve has been shown 
to reproduce the dust reddening of most quasars at z<2.2 
(iRichar ds et al. 2003), while the effective reddening of galaxy 
spectra in the local Univer se is generally mode led by using the 
attenuation law derived bv' Calzetti et al.l (Il994l) . The possibility 
that dust properties may evolve with redshift was first investi- 
gated by Maiolino et al. (2004) who analyzed the spectrum of 
the quasar SDSSJ1048-I-4637 at z=6.2, finding that the redden- 
ing in this object is described better by an SN-type extinction 
curve dTodini & FerraralbOOl b than through the usually adopted 
SMC. A more extended study has been developed by Gallerani 
et al. (2010), who analyzed the optical-near infrared spectra of 
33 quasars with redshifts 4.0 < z < 6.4, finding that all of the 
reddened quasars require an extinction curve deviating from that 
of the SMC, with a tendency to flatten at Aresi < 2000 A. By per- 
forming a simultaneous fit of all reddened quasars in the sam- 
ple, Gallerani et al. (2010) obtained a mean extinction curve 
(hereafter called MEC) at z>4, and the corresponding attenua- 
tion curve (hereafter called MECa,,). The change in extinction 
curve at z>4 suggests that the dust properties at early epochs are 
different with respect to dust at intermediate and low redshifts, 
either because of a different dust production mechanism (e.g. 
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dominated by SNe instead of AGB stars) or because of different 
dust processing (e.g. stronger shocks destroying small grains). 

High-z GRBs provide additional opportunities to study the 
evolution of dust properties at high redshift (e.g. [Stratta et al. 
l2007t lLietal.ll200a IZafar et al.ll2010t iPerlev et al.ll2009.,201Q) . 
In principle, this kind of studies may provide more stringent 
constraints on dust extinction curves relative to QSOs, since 
the intrinsic slope of their continuum can be inferred from the 
associated X-ray data, which are not affected by dust absorp- 
tion. Moreover, GRBs have been discovered at redshifts greater 
than the on es of the most dis tant quasars, e.g. GRB080913 
at z~ 6.7 dG reiner et al. 2009), and GRB090423 at z~ 8.1 
(ISalvaterraet al., 2009: Tanviret al. 2009) . Since (l ong) GRBs 
result from the death of (early) massive stars (e.g. lAbel et al.l 
l2002t ISchneider et al.ll20^ IWooslev & Bloomll2006ir they are 
expected to be detected at even earlier epochs. For these reasons 
GRBs are perfect ca ndidates for investigating both the cosmic 
reionization process (iGallerani et al.ll2008tlPatel et al.l201 0) and 
the evolution of the ISM properties through the cosmic epochs. 

GRB 050904 was a very bright burst at z=6.29 for which 
extensiv e multi-band af terglow follow-up has been performed 



( Cummi ngs et alJ 12005 
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Tagliaferri et al. 2005; Hais lip et alj 



I2OO6: Priceet alJi2006l: fKawai et al..,2006; .Boer et al.ii2o6 6l). At 



the time of its discovery, GRB 050904 was the most distant 
GRB ever detected, and several works have been dedicated to 
broad band data analysis and modeling and to the study of 
the interstellar matter in its host galaxy (e.g. Camp ana et alJ 
120071; iGendre et alll2007h . Dust extinction in GRB 050904 has 
been investigated by several authors, but with discordant re- 
sults. By studying only the spectral energy distribution (SED) 
of the optical/near-IR data, assuming "standard" extinction 
curves (i.e. those of our Milky Way -MW- and for the two 
Magellanic Clouds -SMC,LMC-), several a uthors have found 
very low or zero values of dust extin ction dKann et al. I l2007t 
iHaislip et al.ll2006l: iPerlev et al.ll20To '). From wide-band spec- 
tral afterglow m odeling, ve ry low dust extinction values were 
inferred b y Gou et alj ( |2007|) . assuming an MW extinction curve. 
iLiang & Lii (.200^) . in contrast, find evidence of non negligible 
dust extinction and a best-fit to the data assuming the extinction 
law computed by iLi et aP (2008) with evidence of the 2175 A 
absorption feature. By exploitin g that X-ray em ission is not af- 
fected by dust absorption, Stratt a etakl (l2007h have found evi- 
de nce of dust extincti on at a level comparable to what was found 
by|Liang&Lil(|20Q9!). This result was obtained both through the 
simultaneous optical to X-ray afterglow SED fitting at T-i-0.5 
day and also from the optical/near-IR SED fitting at later epochs. 
These authors also show that an extinction curve expected from 
dust prod uced by supernqvae, s uch as the one found for a z~ 6 
QSO (e.g. lMaiolino et al.ll2004l) . provides a better fit to the data 
than the MW, L MC, SMC extinctio n curves and the Calzetti at- 
tenuation curves, ^afar et al] (l201 d) have published a new care- 
ful data reduction of the photometric data in the Z band, which 
is one of the photometric points most sensitive to dust extinc- 
tion an d to different extinction curves. In their work. lZafar et alj 
(1201 Ol) find that, at Th-0.47 days, the z-band suppression is 
weaker than reported in the previous works and that in gen- 
eral the photometry calibration in this band is affected by sys- 
tematic large uncertainties. By fitting the zYJHKs SED by as- 
suming a power-law spectral model, leaving both the spectral 
index and normalization free to vary, and using any dust extinc- 
tion curve (SMC and the same SN-origin extinction model used 
by Stratta et al. 2007), Zafar et al. (2010) find only marginal ev- 
idence of dust extinction at T-i-0.47 days and no evidence at later 
epochs. 



Given the importance of any confirmation of dust at such 
high redshifts, in this work we carefully revisit the global SED 
analysis for this burst. In particular, contrary to the majority of 
past works based on optical/near-IR data SED fitting alone, we 
exploit the X-ray afterglow emission of this burst to determine 
the spectral slope and normalization of the intrinsic, unabsorbed 
continuum. The difficulty with this burst is that the early X- 
ray afterglow is dominated by strong flare activity. X-ray flares 
are thought to be a distinct component than the afterglow one, 
pos sibly hnked to t he initial burst production mechanism (e.g. 
Falc one et al.ll200'7h . We thus perform a new, detailed, and care- 
ful analysis of the X-ray afterglow of GRB 050904, devoting 
effort to estimating the X-ray flare afterglow contamination at 
the epochs of our interest. We take optical/near-IR data from the 
literature at T-hO.47, 1.25, and 3.4 days, where the photometric 
data in the Z b and is taken from the accurate reanalysis recently 
pubHshed by Zafa r et alj (|2010|) . 

The results of our work are reported in §4 where we pro- 
vide evidence that, at all three analyzed epochs, the optical and 
X-ray afterglow data of GRB 050904 are consistent with being 
extincted by dust, independent of the assumed extinction law 
model. 



2. Swift/XRT data analysis 

Swift/XRT data of GRB 050904 have been calibrated, fil- 
tered, and screened using the XRTDAS package included in the 
HEAsoft distribution (v6.9), as described in the XRT Software 
User's Guid^B- The 0.3-10 keV light curve of this burst is flare- 
dominated at almost all epochs, at leas t up to one day after th e 
trigger time T at Sept 5.0785, 2005 (^Cummings et al.' '2005^ 
Given the flare time-varying energy spectrum (e.g. Falcone et al] 
'2007), time-resolved spectral information was taken into ac- 
count when converting counte-rates to energy fluxes. This proce- 
dure was performed with two methods. In the first one, we used 
the "xrtgrblc" task. This task extracts the light curve from all the 
XRT observations (sequences) performed for this burst (Table 
[TJ, optimizing the source and background extraction regions and 
taking pile-up into account in a time-dependent manner. The 
conversion of the counte-rates into flux units is performed on the 
basis of an energy spectra extraction in temporal bin with width 
that depends from the counte-rate, and on the best-fit of each ex- 
tracted spectrum assuming a double absorbed power-law spec- 
tral model, with one absorption component fixed at the Galactic 
value. Final fluxes are corrected for total absorption. In the sec- 
ond method, rather than extracting the spectrum in a count-rate 
dependent temporal region and fitting the power-law model to 
the data (as the "xrtgrblc" task does), the count-rate to flux con- 
version is performed by assuming the power-law photon index 
that best reproduce the observed hardness ratio and its evolution. 
This met hod is app lied in the "burst analyzer" tooQ (see more 
details in lEvans et al. 2007, 2010). Comparing the light curves 
obtained with the two methods described above, we find over- 
all agreement in terms of unabsorbed fluxes and spectral slopes, 
within errors (Fig[T]). 

However, at late times (t>T+2 days) the uncertainties on 
the estimated flux are quite large. A possible reason may be at- 
tributed to the very low flux level of the afterglow at this epoch: 
while the automatic time-dependent spectral analysis obtained 
with the "xrtgrblc" task and through the burst analyzer tool are 
extremely helpful for the first, bright part of the light curve, they 



http://swift.gsfc.nasa.g0v/docs/swift/a nalysis/l 
|http://www.swift.ac.uk/burst_analyser/| 
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Fig. 1. GRB 050 904 X-ray afterglo w unabsorbed flux as obtained from "xrtgrblc" task (open black circle) and from the "burst 
analyzer" tool by lEvans et alj (1201 Oh (green triangles). The black solid square at about 3.254 days after the burst trigger is the flux 
estimated with the stacking method discussed in the text and used in the temporal analysis. The solid line is the best-fit smoothed 
broken power-law model from the optical data from Tagliaferri et al. (2005) (see inserted panel). The dashed vertical lines indicate 
the three epochs at which the largest spectral optical coverage is available from the literature, namely T+0.47, T+1.25, and T+3.4 
days. The blue dot-short-dashed and red dot-long-dashed lines indicate the expected average flare intensity evolution with time, for 
two different decay rates predicted for late time flares (see §2.1 for more details). The bottom panel shows the X-ray photon index 
as a function of time obtained by both the "xrtgrblc" task (black circle) and the "burst analyzer" tool (green triangles), as well as 
the estimate obtained with the stacking method at about 3.4 days after the burst (black squares). Horizontal dashed line indicates 
the expected value at the epochs of our interest where Vc < vx, for a synchrotron emission from an electron population with energy 
spectral index of p ~ 2.1 - 2.5, that is F^*'' - 2.15 ± 0.10. 



may be not suitable for such faint detection levels. To accurately 
estimate the flux at late times we perform an alternative data 
analysis focused on this part of the light curve. We progressively 
stack the observations starting from T + 2.535 days (sequence 
number 00153514001, see Table [B up to the end of the XRT 
observation at T+9.922 days. The stacking of the cleaned event 
files was performed through the Xselect v2.4a software pack- 
age, while the exposure maps were summed through the Ximage 
v4.5.1 software package. Since the source flux is fading, we 
find a maximum signal-to-noise ratio if we stop our stacking at 
r -I- 4.16 days (sequence number 00153514004). Using the tem- 



poral range found above, we extract the source and background 
energy spectrum from a circular region with a 10 pixel radius 
and an annulus region with an inner radius of 15 pixels and an 
outer radius of 30 pixels, respectively. By grouping the spectrum 
requiring a minimum of 5 counts per energy bin, we obtain a 
source energy spectrum with three energy bins. The model fitted 
to the obtained spectrum is an absorbed power-law model, with 
an absorption component fixed at the Galactic value in the di- 
rection of this burst (4.6 x 10^"cm"^, Kalberla et al. 2005). Past 
results agr ee with no evidence of additional absorp tion at this 
epoch (e.g. lCampana et al.ll2007l: iGendre et al.ll2007h . We find a 
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Table 1. Swift/XRT observations of GRB 050904. From left to right there are the Swift/XRT sequence (observation) identifying 
number, the start and end date of each sequence, and their corresponding distance in days from the burst trigger time. There are no 
XRT observations between 0.6904 and 2.535 days after the trigger. 



Swift/XRT 


T start 


T I 






sequence ID 


YYYY-MM-DD UT 


YYYY-MM-DD UT 


day 


day 


00153514000 


2005-09-04 02:01:26 


2005-09-04 18:26:00 


0.006738 


0.6904 


00153514001 


2005-09-06 14:41:59 


2005-09-06 18:38:56 


2.535 


2.699 


00153514003 


2005-09-07 00:29:33 


2005-09-07 23:32:57 


2.943 


3.903 


00153514004 


2005-09-08 00:48:36 


2005-09-08 05:42:51 


3.956 


4.16 


00153514005 


2005-09-09 00:52:08 


2005-09-09 23:25:57 


4.959 


5.899 


00153514006 


2005-09-10 01:05:23 


2005-09-11 06:17:58 


5.968 


7.185 


00153514007 


2005-09-13 01:03:19 


2005-09-13 23:59:57 


8.966 


9.922 



best-fit photon index of T = 2.4 + 0.5 (68% confidence range) 
and a 0.3-10 keV unabsorbed flux of (2.0+,^°) x lO"''* ergs cm^^ 
s"' (68% confidence range). The obtained flux has been sed at 
the photon weighted mean epoch, that is at T -H 3.254 days, with 
a start time of T + 2.535 days and an end time of T +4.160 days. 

The light curve we use in the following timing analysis is 
the one obtained from the burst analyzer tool since it provides 
a better temporal resolution of the spectral properties. For late 
time observations we use the flux estimated with our faint source 
dedicated analysis as explained above. 

2.1. X-ray flare analysis 

To compare the spectral energy distribution (SED) of the op- 
tical afterglow with its X-ray counterpart, we need the X-ray 
flux level at those epochs with the best optical spectr al coverage, 
which is at T+0.47, T+1.25 and T+3.4 days (e.g. IZafaretal] 
I2010I) . At T+1.25 days no XRT data are available (Table 1). At 
T+0.47 days the X-ray light curve is flare-dominated, while the 
simultaneous afterglow optical counterpart does not show evi- 
dence of such intense flaring activity (FigHJ. There is general 
consensus about c onsidering flares to h ave a different origin than 
the afterglow (e.g. lFalcone et al.ll2007h . therefore we need to es- 
timate the X-ray afterglow flux level excluding the flares. Near 
T+3.4 days there are XRT data, centered at T+3.254 days, thus 
requiring only a small temporal extrapolation. 

Before proceeding to the afterglow flux extrapolation at the 
epochs of interest, we need to know whether late-time (i.e. after 
~ T + 0.5 day) XRT data are still contaminated by flares or not. 
We know that the peak-to-continuum flux ratio AF/ F of X-ray 
flares spans from ~ 1 up to 1000 in some extreme cases (typi- 
cally very early single flares): indeed, the median of AF/ F com- 
put ed over tens of GRBs and hundreds of flares (Chincari ni et al.l 
12010 ) is about 4 and at late times is expected to be 1 . Therefore, 
even if flare activity is present at late times for GRB 050904, 
we expect that the XRT data at T + 3.254 days would repre- 
sent the afterglow flux level enhanced at worst by a factor of a 
few. However, in the following we provide some arguments that 
strongly support the idea that flare activity has deeply damped 
or ceased at late times for this burst. 

Past studies have shown that the average X-ray flare peak 
luminosity < L >, corrected for the contribution of the un- 
derlying afterglow power-law component, decays with time as 
< L >oc r", with a in the range 1 .5- 1.8 (Lazzati et al. 2008). 
Based on a larger sample, recentlv 'Marg utti et al.l (12011) have 
found that < L >oc f-2 7±o i for flares detected up to 1000 s after 
the burst onset (rest frame) and < L >oc f-i-2±o.i j^^gj. ^jj^gg 
Since at late time flares may be affected by a low signal-to-noise 



raticQ, it has been argued that their unbiased intensity decay may 
be steeper For this reason, a joint fitting of early and late time 
flares ha s been performed, pro viding a unique decay index of 
a = 1.8 dBernardini et al.ll201 Ih . 

We exploit these results to have an estimate of the expected 
flare intensity at late times for GRB 050904. The X-ray flux 
starting from about T + 0.08 days (~ T + 1000 s rest frame) 
is completely dominated by flares so that the underlying contin- 
uum does not affect the peak flare evolution with time. In Figure 
1 we plot the average flare intensity evolution, with two possi- 
ble decay rates, where the blue line marks a = 1.8 and the red 
one marks a - 1.2. Each po wer-law has been fitte d to the three 
best-fit flare peaks found by iFalcone et alj ( l2007h at late times 
(i.e. with T + tpeak > 1000 s in the rest frame) for this burst. 
The average flare evolution, computed assuming both a - 1.2 
and a - 1.8, can trace the GRB 050904 flare peak intensity rea- 
sonably well up to about T + 0.5 days, with a preference for 
the shallower decay index (in agreement with past findings, see 
Lazzati et al. 2008, Margutti et al. 201 1). If we assume that flare 
activity is still present at later times, the extrapolation overpre- 
dicts the observed fluxes by more than one order of magnitude. 
In addition, ~ r+0.5 day marks a clear hard-to-soft spectral tran- 
sition, as shown in the bottom panel of Figure 1 . Flares ar e spec- 
trally harder than the underlying continuum (e.g. Margut tTet al.l 
1201 ll) : the rather sharp spectral transition simultaneous to the 
flare peak intensity deviations from the average expected values 
is consistent with a scenario where flare activity has significantly 
damped, or ceased, at epochs later than T + 0.5 days. These find- 
ings enable us to consider the X-ray flux measure at T + 3.254 
days with large confidence as representative of the flare unbiased 
afterglow flux level. 

2.2. X-ray light curve modeling and flux extrapolation 

To extrapolate the X-ray afterglow flux level at the epochs we 
are looking at, we need to assume a temporal model. Past broad 
band modelings of this burst have shown that data from radio 
to X-rays (e.g. iFrail et all l2006i: iGou et all l2007h are consis- 
tent with the fireball paradigms. The afterglow emission is com- 
monly interpreted as synchrotron radiation emitted by a p opula- 
tion of relativistic electrons (e.g. Meszaros & Rees 1997). The 
synchrotron cooling frequency (e.g. Sari et al. 1998) of GRB 
050904, turns out to lie at er iergies below the optical range at th e 
epochs of our interest (e.g. iKann et ani2007t iFrail et al]l2006l) . 
in agreement with the absence of any optical spectral variations 
at those epochs (e.g. Tagliaferri et al. 2005). Thus, X-rays and 



^ Indeed. iMargutti et all i20I Ih have shown that the evolution of the 
flare flux with time resembles the average detectability threshold one. 
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optical afterglow fluxes are expected to decay jointly, following 
the same light curve, with no spectral break between the two en- 
ergy domains. Estimated values of the electron spectral index p 
are within the range 2. 1 - 2. 5 (Frail et al. 2006: iGou et alJl2007t 
IZou et alj|2006l : lchandra et al..,2010) . providing an expected X- 
ray photon index of T"'' = ^"^^ + 1 = p/2 + 1 = 2.15 + 0.10 
consistent with our late time measure. 

The optical afterglow light curve of GRB 050904 is well fit- 
ted by a smoothed broken power-law model between 0.1 and 10 
days after the trigger, with best-fit decay indexes ai - 0.72+0. 10 
and a2 = 2.4 ± 0.4 a n d a break time at ti, - 2.6 +1.0 days 
dTagUaferri et alj|2005h . iKann et all (12007 ') quote ai = 0.85 + 
0.08 and similar 02 and fitting optical data from T + 0.3 days. 
The temporal break is achromatic, and therefore it has been in- 
terpreted as evidence for a jet with opening angle ~ l/r(tt,), 
where r(tb) is the ejecta Lorentz factor at the epoch of the tem- 
poral break (e.g. Sari & Piranll 19991) . 

We thus simply normalize the best-fit optical light curve 
model (from Tagliaferri et al. 2007) at the T + 3.25 day unab- 
sorbed flux value (see Fig.l). The X-ray flux extrapolations at 
T + 0.47, r + 1.25 days, and T + 3.4 days, ai-e 0.03, 0.07, and 
0.001 yuJy. We associated an uncertainty of 25% to the former 
two values and of about 50% to the last value, encompassing the 
uncertainty of the late X-ray flux measures. 



3. Optical to X-ray spectral energy distribution 

Starting from the results obtained in the previous section, we 
take from the literature the optical/near-IR data corrected for 
Galactic absorption at T+0.47, 1.25, and 3.4 days, where the 
photometric data in the Z ba nd is taken from th e accurate re- 
analysis recently published bv'Zafar et alJ (1201 Ol) . and we com- 
pare them with the simultaneous X-ray afterglow unabsorbed 
fluxes obtained in the previous section. We then fit a power-law 
spectral model to the data, setting the spectral index and normal- 
ization to the values obtained from X-ray afterglow. No spec- 
tral variation is expected during each epoch since Vc < vx^opt at 
these times (see §2.2). Given the large uncertainty affecting the 
X-ray photon index we have estimated at T + 3.254 days, we 
conservatively performed our analysis using both the expected 
(F'^'' = 2.15 + 0.10, 13'^'' = 1.15 ± 0.10) and the measured 
(Fx = 2.4 + 0.5, /5x = 1.4 + 0.5) X-ray photon indexes. 

Figure |2] shows the obtained broad band SEDs at the three 
epochs. Squares give the X-ray fluxes estimated by us in the 
previous se ction, while diamonds are the values reported in 
IZafar et aP ([2010). We note t hat our X-ray fluxe s deviate signif- 
icantly from those reported in lZafar et al.l (|2010) at T + 0.47 and 
r+ 1 .25 days. The latter discrepancy may be due to different l ight 
curve modeling and flux extrapolation. In lZafar et al.l (1201 Ol) the 
X-ray light curve is modeled with a smoothed broken power- 
law, however parameter values are not reported, hence prevent- 
ing further comparative studies. The hatched areas (delimited by 
dashed lines) indicate the range of intrinsic power-laws consis- 
tent with the +2cr uncertainty on /3x- The dashed line in the mid- 
dle shows the power-law resulting from the best-fit slope to the 
X-ray data at late epochs (j3x = 1 .4). The gray shaded areas show 
the range of intrinsic power-laws consistent with the +lcr uncer- 
tainty on fix- The cyan shaded areas show the power-laws ex- 
pected from synchrotron emission (for an electron spectral index 
~ 2.1 - 2.5 and Vc < vx^opt, see §2.2). The colored areas show 
the optical to X-ray SED best-fit obtained assuming a power-law 
model with index free to vary within the range /3x ± 2cr (Tab. 3) 
and the extinction/attenuation curve associated with each panel. 



To quantify the need of dust extinction and reddening we 
have fitted the data at each epoch with different extinction and 
attenuation curves. The question of whether dust reddening is 
better described by an extinction curve or by an attenuation 
law depends on the geometry of the system. For quasars, the 
simple dusty "screen" geometry applies, but for galaxies one 
has to consider that generally dust is mixed with the emitting 
sources (either stars or ionized gas). As a consequence, atten- 
uation curves i.e. the ratio of the observed spectrum to the in- 
trinsic total light emitted by the whole system before dust ab- 
sorption, are likely more appropriate for galaxies. In the case of 
GRBs it is not clear which geometry is more appropriate. The 
screen case is certainly more appropriate when the size of the 
emitting optical afterglow is smaller than the distribution of the 
dusty medium. However, the dust may well be the same as pro- 
duced by the progenitor before the explosion. In this case the 
size of the dusty medium can be comparable to the one of the 
expanding afterglow, and an attenuation curve may be more ap- 
propriate. As a consequence, to describe dust reddening, we con- 
sider the SMC extinction curve, the Calzetti attenuation lawQ, the 
mean extinction curve (MEC)_resulting from the analysis of 33 
quasars at z> 4 done by iGallerani et alJ (120101) . and the corre- 
sponding attenuation curve ( MECqtt). We also test th e SN-type 
extinction curve proposed bv lTodini & Ferraral (1200 lb . which re- 
produces the dust extinction observed in a BAL QSO at z=6.2 
(Maiolino et al. 2004; Gallerani et al. 2010). In Figure|3] we plot 
the extinction and attenuation curves adopted in our analysis, 
normalized to Asooo, i.e. to the extinction value at the rest frame 
wavelength of Aresi = 3000 A. 

We report the best-fit results for each extinction/attenuation 
curve and at each epoch in Tables |2] |3] and |4] where we give 
the best-fit /5x and A3000, as well as the resulting best-fit re- 
duced and the associated probability. Within each table we 
leave /5x free to vary in different intervals: /3x £ [0.9-1.9] in 
Table |2] and /3x £ [0.4-2.4] in Table [3] which correspond to the 
Icr and 2cr errors on the measured X-ray photon index at late 
times, respectively (see §2.1); fix e [1.05-1.25] in Table|4]is the 
range expected from synchrotron emission with < vx,opt and 
p ~ 2.1 - 2.5 (see §2.2). The data in bold highlight the ex- 
tinction/attenuation curve giving a good fit to the data {P{x^ < 
x\p) < 68%) at aU epochs. 

In some cases the fix best-fit value coincides with the lower 
extreme of the interval considered. In these cases, we cannot 
compute the Icr error on Asooo and we can only provide an in- 
terval within which this parameter is allowed to vary. Moreover, 
for some combinations of fix intervals and extinction curves, the 
best-fit values are characterized by a P(x^ < x\f) > 68%; in 
these cases, we only provide the best-fit values. 

4. Results 

Tables |2] to |4] show that, regardless of the interval adopted for 
fix and taking only the best-fit with P(x^ < x^^) < 68% into 
account, all epochs require dust extinction. The SMC and SN- 
type extinction curves provide a good fit to the data at all the 
three analyzed epochs (Fig. 2). 

At r + 0.47 days and at T + 1 .25 days the obtained extinction 
values are systematically higher and lower, respectively, than 
those obtained at T + 3 .4 days. Since such an evolution with time 
has no physical explanation, we interpret these discrepancies as 

It has been shown that the Calzetti law is the result of the attenu- 
ation by a medium with an SMC-type d ust dOordon et alJI 19971: llnoud 
2005). But see also lPanuzzo et all d2007h and lPierini elail j2004r) 
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Table 2. Best-fit model parameters obtained for l3x,opt - Px ^ [0.9; 1.9], i.e. Icr interval of the X-ray energy spectral index value. 

Epoch Extinction curve Px,opi A3000 Xbf ^(x^ < xIf) 

0.47 days SMC 1.15!^{;* 036^^ 5:92 (UO 

MEC 1.19 0.69 1.68 0.81 

Calzetti 1.22 0.99 1.35 0.74 

MEC„„ 1.29 1.50 2.37 0.91 

SN-type 1.16;!ig^ ^-^^Tn] ^-^^ ^'^^ 

1.25 days SMC 1.11+™ OW^^ OiOT 0:02 

MEC i-i2!oll o.i4!;;;j« 0.09 0.04 

Calzetti 1-13;!! 11 0.22;" 0.08 0.03 

MEC„„ l-15!g|l <0.43 0.11 0.04 

SN-type l.li;^-;] O-H 

3.4 days SMC 1.12+|* 0.19+^-^ 03 02 

MEC l-15!nn! 0-38!o-n7 0.30 0.26 

Calzetti l-17!oi8 0.57!^°^ o.22 0.20 

MEC„„ l-^l^nns 0-90!o28 0.46 0.37 
SN-type i-i^-lm O-^'^-nm 0-0^ "-O^ 

Table 3. Best-fit parameters obtained for /3x.opt - Px ^ [0.4; 2.4], i.e. 2cr interval of the X-ray energy spectral index value. 



Epoch 


Extinction curve 




^3000 


Xbf 


P(x'<xiF) 


0.47 days 


SMC 


1 |s;+0.04 
^■^-"-0.04 


36+" "' 

"■-'"-0.01 


0.92 


0.60 




MEC 


1.19 


0.69 


1.68 


0.81 




Calzetti 


1.22 


0.10 


1.35 


0.74 




MEC„„ 


1.29 


1.50 


2.37 


0.91 




SN-type 




51+""^ 


0.02 


0.02 



1.25 days SMC l.ir^- j Oi)^™ OiOT (0)2" 

MEC l-12!n 1 0-1441 0.09 0.04 

Calzetti l-13!o|i 0.22^™^ 0.08 0.03 

MEC„„ 1-15!!!|2 <0-43 0.11 0.04 

SN-type 1.11+°- 1 ■ 0.09+"";! 0.11 0.04 

3.4 days SMC 1.12!^* Ol95f 0^3 OT 

MEC l-15!nn! 0.38!|| 0.30 0.26 

Calzetti l-17!oi8 0.57!^^^ o.22 0.20 

MEC„„ l-^l^fiftg 0.90!™^ 0.46 0.37 

SN-type 1.13+°-i;^ 0.27+;;-°3 0.04 0.04 



Table 4. Best-fit model parameters obtained for /3x,opi - Px ^ [1 .05; 1 .25] (i.e. the interval of the theoretically expected X-ray photon 
index for an electron spectral index p ~ 2.1 - 2.5 and Vc < vx-) 



Epoch 


Extinction curve 


Px,opl 


^^3000 


Xbf 


P(x'<rBF) 


0.47 days 


SMC 

MEC 


1 |j:+0.04 
^■'--^-0.04 

1.19 


36+" "' 

"--'"-0.01 

0.69 


0.92 

1.68 


0.60 

0.81 




Calzetti 


1.22 


0.99 


1.35 


0.74 




MEC„„ 


1.25 


1.39 


2.60 


0.92 




SN-type 


1 16+""' 


Q g,+0.03 
"--'^-0.03 


0.02 


0.02 


1.25 days 


SMC 


1.11+"'| 


09+" "' 

-MZ 


0.07 


0.02 




MEC 


J j2+8:li 


0.14+"" 


0.09 


0.04 




Calzetti 




22+''-^'^ 


0.08 


0.03 




MEC„„ 
SN-type 


^-^-^-0.12 
1 11+"" 


42+" "2 

"-^^-0.31 

11 +""2 


0.11 
0.11 


0.04 
0.04 


3.4 days 


SMC 

MEC 


^-^^-0.08 
1 1 c+0.08 


0-194! 

0.38!"!^ 
0.57+""^ 


0.13 

0.30 


0.12 

0.26 




Calzetti 


J j-y+O.OSi 


0.22 


0.20 




MEC„„ 


J 21+^-^^ 


0.90+1*1 


0.46 


0.37 




SN-type 


^■^■'-n 118 


27+""° 


0.04 


0.04 



the result of the systematics affecting the optical data reduction 
(see Zafar et al. 2010). By averaging the results obtained assum- 
ing the SMC and the SN-type extinction curve over all epochs, 
the amount of dust absorption at Ai^st - 3000 A is at a level of 
< ^3000 >= 0.25 + 0.07 mag where the uncertainty takes our 
ignorance on the true extinction curve into account. The average 
best-fit optical to X-ray spectral index, independent of its initial 



fixed range (Tables 2, 3 and 4), is < Px.opi >- 1-13 ±0.22, which 
is nicely consistent with the range of values expected from past 
broad band modeling of this burst, which is pjl with the elec- 
tron spectral index in the range p ~ 2.1 - 2.5 ( Frail et al.. ,2006t 
IGouetal.1120071: IZou et aWmM IChandra et al.l2010h . 

We find that the average properties of X-ray flares 
(iMargutti et all l20lU iBernardini et al.l 120111: IChincarini et al.1 
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Fig. 2. Optical to X-ray SED at 0.47, 1.25 and 3.4 days after the burst trigger (from left to right) and fitted using the two extinction 
curves that provide excellent fit at all epochs, that is, the SMC and the SN- type (from top to bottom). Squares show the X-ray flux 
obtained by our analysis while diamonds are the fluxes bv lZafar et al.l ( l2010h . The hatched area (delimited by dashed lines) indicates 
the range of intrinsic power-laws consistent with the +2cr uncertainty on fix- The dashed line in the middle shows the power-law 
resulting from the best-fit slope to the X-ray data at late epochs {fix - 1.4). The gray shaded areas show the range of intrinsic 
power-laws consistent with the +lcr uncertainty on [ix- The cyan shaded areas show the the power-laws expected from synchrotron 
emission (for an electron spectral index p ~ 2.1 - 2.5 and Vc < vx,opt)- The colored areas show the optical to X-ray SED best-fit 
assuming a power-law model with index free to vary within the range Px ± 2cr (TableO. 



1201 0') provide convincing indications that flare activity is 
strongly damped or has ceased at late times for GRB 050904. 
We estimate the transition from flare-dominated to afterglow- 
dominated X-ray flux at about T + 0.5 day, where the expected 
average flare peak intensity starts to overpredicts the observed 
fluxes and when an abrupt hard-to-soft spectral transition is ob- 
served (Fig. 1). The X-ray flux measure at T-i-3.254 days is more 
than one order of magnitude lower than the expected flare inten- 
sity, therefore likely representative of the flare-unbiased after- 
glow flux level. Results quoted in Tables 2,3 and 4 were obtained 
by extrapolating the X-ray afterglow at the three epochs where 
best optical/near-IR spectral coverage is available and by nor- 
mahzing the best-fit temporal m odel obtained at optical wave- 
lengths ('Tag liaferri et al.l r2005h at the X-ray unabsorbed flux 
value measured at T H- 3.254 days (Fig.l). 

It may be argued that several GRBs show that X-ray and 
optical after glow hght curves do not decay with the same be- 
havior (e.g. iPanaitescu et"ani2006l) . However, even though there 



is also a non negligible fraction of GRBs for which the opti- 
cal and X-ray emission decay jointly (e.g. Gates et al. 201l|, 
GRB 050904 is one of the few bursts for which a broad band 
modeling was feasible from radio to X-rays. In particular, it 
has been found that data are generally consistent with fireball 
paradigms, with the synchrotron cooling frequency below the 
optical range at the epochs of interest, thus supporting our as- 
sumption (e.g. lFrail et al.ll2006l) . In addition, results for dust ex- 
tinction at r H- 0.47 and T + 1.25 days are on average consis- 
tent with the A3000 values obtained ?A T + 3.4 days, when the 
S wift/XRT observations are available at T H- 3 .254 days, and thus 
flux extrapolation depends little on the assumed temporal model. 
At r -H 3.4 days, the extracted X-ray to optical SED shows clear 
evidence of dust extinction, with A3000 in the range that goes 
from 0.2 to 0.9 mag, depending on the assumed dust recipe (Tab. 
4). 

These results indicate, beyond a doubt, that the primeval 
galaxy at z=6.3 hosting this GRB has already enriched its ISM 
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Fig. 3. Extinction and attenuation curves adopted in our analysis, 
normalized to the extinction value at the rest frame wavelength 
ofArest = 3000 A. 

with dust. However, while the presence of dust attenuating the 
GRB 050904 optical afterglow at any epoch is firmly estab- 
lished, the type of extinction/attenuation curve is not well con- 
strained, although we find that only the SMC and SN-type ex- 
tinction curves provide a good fit to the data {P{x^ < x\f^ < 
68%) at all epochs (Fig. 2). 

5. Summary and conclusions 

We reanalyzed the afterglow of GRB 050904 at z=6.3 at those 
epochs where the best spectral coverage is available in the 
optical/near-IR range (UV rest frame), namely at 0.47, 1.25, 
and 3.4 days after the trigger, by fitting the simultaneous 
optical/near-IR and X-ray SEP. In thi s work , we exploited the 
recent reanalysis done by IZafar et aTl (1201 Oh on Z-band data, 
which are the most sensitive to dust reddening for this burst. 

We analyzed the Swift/XRT data of GRB 050904, taking 
special care to estimate the flare unbiased afterglow flux level 
at the epochs of interest and to extrapolate the flux where no 
data are available. A spectral power-law model was fitted to 
the extracted SEDs, setting the spectral index and normaliza- 
tion at the dust-unbiased X-ray emission values. No spectral 
breaks or chromatic temporal breaks between X-rays and op- 
tical wavelengths are expected at these epochs, according to past 
broad band mo deling results for this burst (e.g. lFrail et al.ll2006t 
iGou et af]|2007l) . Indeed, this predicts that the synchrotron cool- 
ing frequency Vc is already below the optical range at the time 
of interest, thus an identical afterglow decay law at X-ray and 
optical wavelengths is expected. 

We investigated any presence of dust extinction in the GRB 
afterglow by using the SMC extinction curve and the Calzetti 
att enuation law. We also used the MEC at 4.0<z<6.4 inferred 
by iGallerani et al.l (1201 Oh from the analysis of 33 quasars, and 
the associated attenuation curve (MEC,,,/). The SN-type extinc- 
tion curve, proposed by Todini & Ferrara (2001), which repro- 
duces the dust extinction observed in a BAL QSO at z=6.2 
(iMaioUno et al.ll2004t iGaUerani et al.ll2010l) was also tested. 

From a simultaneous fit of the rest frame UV to X-ray SED, 
we find clear evidence of dust absorption at all epochs. The SMC 
and SN-type extinction curves provide good fits at all epochs, 
with an average value of dust extinction of 0.25 + 0.07 mag at 
Arest - 3000 A. At r -I- 3.4 days where nearly simultaneous X- 



ray and optical data are available, thus where X-ray temporal ex- 
trapolation weakly depends on the assumed model, the extracted 
X-ray to optical SED shows clear evidence of dust extinction, 
with A3000 in the range that goes from 0.2 to 0.9 mag, depending 
on the assumed dust recipe (Tab. 4). 

Our findings indicate that the primeval star-forming galaxy at 
z - 6.3 hosting this GRB has already enriched its ISM with dust. 
The type of extinction/attenuation curve is not well constrained 
by the data, although we find that only the SMC and SN-type 
extinction curves provide good fit to the data at all epochs. 

We emphasize that our results were obtained from a simul- 
taneous optical to X-ray SED fitting, while most of the previous 
studies were performed on the optical/near-IR SED alone. Our 
method strongly benefits from the information from past broad 
band modeling and from the X-ray afterglow flux normalization 
level and spectral slope at 3.25 days after the trigger, when we 
demonstrated that the flare activity has very likely ceased. 
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